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Abstract 

Small momentum expansion of the "sunset" diagram with three different masses 
is obtained. Coefficients at powers of p 2 are evaluated explicitly in terms of diloga- 
rithms and elementary functions. Also some power expansions of "sunset" diagram 
in terms of different sets of variables are given. 



The "sunset" diagram (see Fig. 1) was the object of investigation in several recent 
works 0, H H §], [5]. Such activity is initiated, certainly, by rather high precision of 
experimental testing of the standard model that needs sometimes evaluation of two-loop 
Feynman diagrams for comparison of the theory and experiments. 




Figure 1: The "sunset" diagram 
*e-mail address: lunev@hep.phys.msu.su 



1 



Probably, the "sunset" diagram cannot be expressed through known special functions 
of one variable. So for investigation of this diagram it was used numerical methods J|. H 
and some approximations schemes, such as expansions in powers of p 2 (and possibly 



In particular, in paper || it was obtained small momentum expansion of "sunset" 
diagram in terms of Lauricella series. These series are convergent in the region 



where m^, i = 1,2,3, are masses in propogators. But for majority of applications 
mi + m2 > m-3 for any numbering of masses. Such situation takes place, for instance, in 
the course of evaluation of self-energy of Higgs boson (when mi = 777,2 = 777,3 = mil), and 
in the course of evaluation of contribution to Z-boson self-energy due to vertex Z 2 WW* 
(when mi = 7772 = raw, "13 = m^). 

The aim of the present paper is to complete the results of and to obtain small 
momentum expansion of "sunset " diagram that can be used for all values of masses, and 
to evaluate in closed form coefficients at powers of p 2 . In addition, we will give some 
power expansions of "sunset" diagram that are convergent, at least, in the region 



Obviously, these expansions can be used for all values of masses if one defines m 3 as 
min{mi, 7722, 777,3}. 

Let / = I(p 2 , mi, m 2' m V) b e ^ ne "sunset" diagram. Our main tool will be the following 

one-fold integral representation for the function J: 

1 See also an old paper by Mendels [|| where it was obtained expansion of "sunset" diagram in equal 
mass case in powers of p 2 /p 2 + m 2 
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where J is one-loop "bubble" diagram (see Fig. 2), p is spectral function of J, 



P = TV 
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and functions /j, z = 1, 2, 3 must be defined in such a way that integral (^j converges. 



(4) 





Figure 2: The "bubble" diagram 

Representation (||P was obtained in a previous author's work 0. Very similar integral 
representation was obtained independently in above mentioned papers [|3], f|]. Another one- 
fold integral representation was given in .An analogous one-fold integral representations 
for the five propogator self-energy diagram was derived in || HI. 

One notes that, due to renormalization freedom, function I is defined up to polynomial 
of the first degree with respect to p 2 . So it is sufficient to evaluate the function 



dl 



Comparing 



I(p 2 ,...) = I(p 2 ,...)-I(Q,...)-p 2 ^(0,...) (5) 

ap z 

and (El), one obtains the following integral representation for /: 



2 In |Q] representation (|^) was derived only for the case mi = TO2, but, in fact, this condition was never 
used and for general case representation (0) can be obtained in the same way. 
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where 



(6) 



dJ 



J(p\...) = J(p 2 ,...)- J(0,...)-P 2 ^(0,...) 
The function J can be represented by the following power series: 



(7) 



J=7r2 (k+i-i)\(k+i)\ ( P 2 y ( x m^' 
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The proof of @ is very simple One considers the representation of J through 
Feynman parameters: 



J=-n 2 /^log 
Jo 



p 2 gi - Q + frnj + (1 - j)a 2 
a 2 



Expanding the integrand in powers series with respect to 



V\{1 - Q + Ornl - a 2 ) 
and integrating, one obtains (H). These operations are correct if 



max 

0<£<l 



(T- 



< 1 



From (|K]) one can derive that series @ converges, at least, if 



(9) 



(10) 
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|p | + m 3 < cr 
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One introduces dimensionless variables 



3 Another proof was given in p 
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(mi + m, 2 ) 2 mj 

a 2 " (mi + m 2 ) 2 ' 



(mi + m2) 2 ' (mi + m 2 ) 2 



In these variables, using definition of Gauss hypergeometric function 2-^1, one can 
rewrite (|j) in the following way: 



J = tt 2 Y: x k z k \ \ ; 2 F x {k, + 2fc + 2; 1- xy) (13) 

But 



(fc-l)!fc! , , s 1 d k , sh . , d k log(l - u) , , 

(see, for instance, pTT[]) . From the last formula, after some elementary transformations, 
one can obtain: 



\ 2k + 1)l ^(k, k + 1; 2k + 2; 1 - xy) 
1 1 d d k k+1 d k ~ l \og(xy) 



k\(k + l)\x k 2 dx dy kU dy k 1 y(l — xy) ^ 
Substituting (|13"D and (|T5|) in (||) and integrating by parts with respect to x, one finds: 



z k 



where functions fk(y, A) are defined by formulas: 



tt ^ d k k+1 d*' 1 f(y, A) 
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/(y,A)=logy+ f^^(l-,)(l-A|/) (18; 
jo 1 — xy ax y 

The integral in flT8|) can be easily evaluated and we obtain for f(y,X) the following 
explicit expression: 



/(y,A) = logy- 



Ay 



-Sp 
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where Sp(x) is the Spence function (or dilogarithm) , 



Sp 



1-t 2 



(19) 



i + Va 

1-vT C y(l-A) 



(2A - y(A + 1) + 2 v /A(l-y)(A-y)) 



(20) 



Formulas (Tl6|), (|iT|), and (|19|) give explicit expressions for coefficients at powers of 
p 2 in Taylor expansion of "sunset" diagram. But formula flTQ] ) seems too complicated. 
Fortunately, in almost all applications ui\ = m 2 (that is, A = 0), and just in this case 
formula (|19|) considerably simplifies: 



f{y,0) = log y 



where Chi^O) is Clausen integral. 



y/vO- - v) 



Cl 2 (2 arcsin y/y) 



(21) 
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In conclusion, we will derive two series expansions for "sunset" diagram for the case 
of three different masses, when direct using of explicit formulas (0), (|17D, and ([19]) is 
difficult due to complexity of formula ([19|) . 

First, one substitutes the decomposition (H) in (H), expands integrand with respect to 
A and performs integration. Then, after some elementary transformations, one obtains 
the following result: 



ttV) 2 / fk + l + n + l)\(k + l + n + 2)\T(n + l) 



mi + m2 ) 2 i,M>o \n?o n\(2k + l + n + 5)\T(j + k + l + n + 



x (^g£Xpi)v(i- s) v (22) 



Above mentioned operations are correct if condition flT5|) is valid for all a > (mi + 



m2) 2 - This means, that series (|22|) converges, at least, at region @. 



The sum over n can be evaluated in terms of finite sums. But the corresponding 
expression is too cumbersome to be useful. So, if it is desirably to have closed expression 
for coefficients at X^y l z h , it is preferable to use another series expansion for I, namely: 



t 4, _l v- **A* lj (j + k-2)\ 



j>0,fe>2 



2j! u 2 ; 1 k(k + l)T(j + k + ±] 



(k + l)\(k + l + l)\(j + k + l-l)\ 



where 



hjki = ^{j + k + I) + ip(k + I + 1) + ip(k + I + 2) 



-i/>{j + k + l + ^)-iP{l + l)-il>(l + 2) 
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In order to prove (^), it is sufficient to write function 2-F1 in ([13]) as power series with 
respect to y and log y, to substitute this series in (|6|), and to perform the integration. 

These operations again are correct, if condition (11) is valid for all a 2 > (mi + m.2) 2 . 
So the region of convergence of series fl23f) is not less than one defines by formula @). 

Large momentum expansion for "sunset" diagram also may be derived from integral 
representation (|3|). But in this case results of seem quite exhaustive and so further 
investigations of this problem are unnecessarily. 

Author is indebted to A.I.Davydychev and B.A.Smirnov for interesting discussions 
and comments. 
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